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PREFACE 


The  author  of  this  report  la  grateful  to  the  following  project 
engineer*  at  RADC  (IRAP),  Griffis*  AFB  NY  1 3441  for  their  advice  and 
aaalatance  In  generating  this  report:  David  Clark,  Albert  Proctor, 

Walter  Csytyckl,  and  Bruno  Beek. 

Several  government  agencies  (Naval  Electronic  Laboratory  Center) 
and  private  corporations  (American  Electronic  Laboratories)  are  known  to 
be  using  Walsh  Transform  processing.  However,  up  until  now,  very  little 
information  has  been  published  describing  the  technology  Involved. 
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SECTION  I 


Introduction 

In  the  area  of  digital  signal  processing,  the  representation  of  Cine 
domain  data  by  the  superposition  of  simple  functions  Is  an  Important  analysis 
technique.  This  Is  often  done  in  the  frequency  domain  through  the  use  of  the 
Fourier  transform.  Once  the  transform  of  the  input  data  13  obtained,  the 
Fourier  coefficients  which  contribute  to  the  Information  part  of  the  signal 
can  be  retained,  while  those  that  contribute  only  to  the  noise  can  be 
dropped. 

In  the  Fourier  transform,  sine  and  cosine  terms  for  the  set  of  basis 
functions  which  are  used  to  selectively  reconstruct  the  Information  contained 
in  the  Input  data.  Most  naturally  occurring  disturbances  are  sinusoidal  In 
nature,  example:  all  forma  of  electromagnetic  radiation,  human  speech, 
seismic  events,  etc.,  and  thus  the  Fourier  transform  has  proven  to  be  useful 
In  the  analysis  of  these  signale.  Many  forms  of  Information  collection  and 
transmission  systems  built  by  man,  however,  do  not  consist  of  sinusoidally 
varying  values,  but  of  rectangularly  shaped  on-off  pulses.  Some  exa spies  are 
digital  coanunlcatlons  channels,  radar  pulse  envelopes,  and  telemetry  data. 
The  analysis  of  these  signals  through  the  utilisation  of  a transform  that  Is 
based  on  sine  and  cosine  functions  can  be  cumbersome. 

Walsh  functions  form  s complete  orthonormal  set  of  rectsngulsr  waves. 

The  first  eight  ere  >hown  in  figure  one.*  There  are  several  ways  to  generate 
(mathematically)  a Walsh  function  series.  Each  has  certain  advantages: 

1.  By  means  of  s difference  equation 

2.  From  products  of  Rademacher  functions 

3.  Through  the  Hadamard  matrix 

4.  By  the  use  of  Boolean  Synthesis. 


the  derivation  of  Walsh  functions  Is  not  the  subject  of  this  study  and  more 
information  In  this  area  can  be  obtained  from  reference  1 at  the  end  of  this 
report. 

Walsh  function*'  take  on  two  possible  values,  +1  or  -1.  The  Walsh 
transform  can  be  used  to  process  time  domain  signals  of  any  form,  but  are 
particularly  useful  for  signals  which  are  rectangularly  shaped. 
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SECTION  II 

Comparison  of  Walsh  and  Fourier  Spectrum  Processing 

The  decision  of  whether  to  use  Fourier  or  Walsh  domain  processing 
techniques  Is  dependent  upon  the  characteristics  of  the  signal.  A smoothly 
varying  continuous  signal  favors  the  Fourier  transform  whereas  a rectangular 
waveform  Is  reconstructed  with  fewer  terms  by  using  the  Walsh  transform. 

(See  figure  2) 

The  Fourier  transform  exhibits  the  Gibb's  phenomena  (divergence  at  the 
aero  crossing)  and  this  necessitates  the  use  of  a large  number  of  high  frequen- 
cy terms  in  the  Fourier  transformation  to  retain  the  square  corner 
characteristics  of  a pulse.  The  Walsh  transform  is  uniquely  adapted  to  pulsed 
signals  and  does  not  exhibit  the  Gibb's  phenomena. 

The  Walsh  transformation  Is  a real  one,  while  the  Fourier  transformation 
la  complex.  In  digital  signal  processing,  the  Fast  Fourier  Transform  (FFT) 
requires  a total  of  N LOG^N  complex  multiplications  and  additions.  The  Fast 
Walsh  Transform  (FWT)  requires  N LOG^N  real  additions  or  subtraction*.  This 
makes  the  FWT  much  faster  than  the  FFT.  Also,  sine  and  cosine  functions 
cannot  be  represented  exactly  by  a finite  number  of  bits;  thus  a source  of 
truncation  noise  Is  Introduced  by  the  discrete  Fourier  transform  which 
involves  repeated  multiplication  by  a complex  number.  The  FWT  Involves  only 
additions  and  subtractions  of  real  numbers  and  precise 
possible.  Thus,  the  FWT  Is  not  a source  of  noise. 

The  amplitude  coefficients  of  the  discrete  Fourier 
affected  by  the  phasa  of  the  Input  signal.  This  causes 

the  same  complex  amplitude  spectral  coefficients  regardless  of  its  position 


representation  are 


transform  are'  not 
a given  signal  to  have 
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FIGURE  2 Synthesis  of  Pulsed  Signals 
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within  the  sampling  window.  The  Walsh  transform  will  exhibit  a different 
spectrum  as  the  signal  is  shifted  within  the  window.  (See  figure  3) 


SECTION  III 


Data  1 Description 

In  each  of  the  data  plots  that  follow  a 256  point  window  Is  used. 

Present  in  this  window  are  one  or  acre  pulsed  signals  plus  added  Caussian 
noise.  Pulses  of  different  widths  and  slgnal-to-noise  ratios  are  examined. 

Data  Table  1 lists  all  the  parameters  which  pertain  to  data  plots  1-1 
through  1-13.  Only  low  pass  filtering  of  the  Walsh  spectrum  is  performed. 

Line  1 of  each  data  page  shows  the  input  time  domain  pulse  plus  added 
Caussian  noise. 

Line  2 is  the  Walsh  transform  spectrum  of  the  line  1 input  data. 

Line  3 shows  the  coefficients  which  were  used  eo  reconstruct  the  time 
domain  input  data.  This  is,  in  effect,  sequency  domain  low  pass  filtering. 
Only  the  terms  that  are  plotted  in  this  line  are  used  in  the  reconstruction; 
all  others  are  set  equal  to  aero. 

Line  4 is  the  reconstructed  time  domain  data.  It  was  obtained  by  taking 
the  Inverse  FWT  of  line  3. 

Line  5 is  the  time  domain  amplitude  threshold  of  line  4.  All  data 
below  the  given  threshold  value  are  not  output  or  plotted. 

Clarification  of  the  Data  Table  1 entries  is  given: 

START  OF  PULSE  - This  gives  the  number  of  the  sample  in  the  256  point 
window  where  the  pulse  begins. 

END  OF  PULSE  - This  is  the  number  Of  the  sample  where  the  pulse  ends. 

This  information  is  given  because  it  will  not  always  be  obvious 
in  low  signal  to  noise  ratio  data.  The  pulse  position  is  also 
outlined  by  dashed  lines. 
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SEQl'ENCY  LOW  PASS  it  OF  TERMS  - The  Walsh  coefficients  are  plotted  In 
ascending  sequency  order.  The  number  of  low  order  terms  used 
when  the  Inverse  FWT  Is  calculated  Is  given  here.  All  other 
coefficients  In  the  array  are  set  equal  to  zero.  This  Is  the 
line  3 data. 

INPUT  DATA  S/N  (db)  - The  signal  to  noise  ratio  of  the  line  1 Input 

data  is  given  here.  There  are  many  methods  of  calculating  S/N 
ratios.  Here  the  Input  signal  Is  assumed  to  be  a voltage  and 
the  formula  used  is: 

(average  value  of  signal) 

S/N(db)“  20  LOG  (average  value  of  noise) 

TIME  THRESHOLD  - The  last  line  (line  S)  of  the  data  plota  shows  how  all 
of  the  reconstructed  data  below  the  listed  value  of  amplitude 
threshold  in  the  time  domain  la  truncated.  Only  the  samples 
which  have  values  greater  than  this  threshold  value  are  plotted 
on  the  screen.  This  cleans  up  the  pulse  further. 

, t 

' i 

PULSE  WIDTH  ERROR  LEADING  EDGE  - A low  S/N  ratio  pulse  will  often  | 
exhibit  some  error  as  to  the  exact  point  where  the  pulse 
begins.  The  number  of  samples  that  the  leading  edge  of  the 
reconstructed  pulse  is  In  error  as  compared  to  Its  position 
in  the  line  1 Input  data  Is  given  here. 

PULSE  WIDTH  ERROR  TRAILING  EDCE  - The  number  of  samples  that  the 
.reconstructed  pulse's  trailing  edge  is  in  error  as  compared 
to  line  1 Is  given  here. 


9 


DATA  TABLE  1 


Notes  on  Data  Table  1 


Data  1-1  and  1-2  - These  plots  have  identical  input  data,  but  1-1 
processes  30  Walsh  coefficients  In  order  to  obtain  the  reconstructed  tne 
domain  pulse  while  1-2  processes  only  8.  In  this  case,  8 tarns  reconstruct 
the  pulse  with  the  same  accuracy  of  pulse  width  and  amplitude  as  do  30  terms. 
The  first  4 sequency  terms  dominate  because  of  their  large  amplitudes. 

Data  1-4  - Two  pulses  of  differing  pulse  widths  exist  in  the  window. 

The  first  pulse  has  a width  of  30  samples.  The  second  pulse  has  a width  of 
IS  samples. 

Data  1-11  - At  this  low  S/N  ratio  (-4  db)  only  one  pulse  actually  exists 
in  the  noise.  Its  position  is  indicated  by  the  dashed  lines.  The  other 
pulse  that  shows  up  in  line  5 is  a false  indication.  A close  look  at  the 
input  data  reveals  that  the  noise  does  take  on  slightly  more  positive  values 
where  the  pulse  is  indicated. 

Dat*  1-12  and  1-13  — These  plots  have  identical  line  1 input  data.  It 
consists  of  3 pulses,  10  samples  wide,  with  10  samples  separation.  Data 

line  5 Incorrectly  shows  only  one  pulse  present  in  tha  noise.  Data 
1-13  shows  all  3 pulses  because  a sufficient  number  of  terms  was  used  in 
the  reconstruction  of  the  input  dsts.  If  sn  insufficient  number  of  term  is 
processed,  several  closely  spaced  pulses  nay  be  mlslntrepreted  and  combined 
into  one  pulse,  particularly  at  low  S/N  ratios. 
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SECTION  IV 


Data  2 Description 

Data  Table  2 llata  all  the  parameters  which  pertain  to  data  plots 
2-1  through  2-7.  In  this  data  set,  no  low  peas  filtering  of  the  Welsh 
spectra*  was  performed.  Instead,  amplitude  thresholding  of  the  entire 
Walsh  spectrum  takes  place.  All  sequencer  terms  which  have  amplitude  values 
below  the  listed  threshold  are  set  equal  to  sero.  This  threshold  value 
was  vsrled  from  .3  to  1. 

Notes  on  Data  2 


Data  2-1,  2-2,  and  2-3  - These  data  plots  have  Identical  Input  data. 

In  each  one,  however,  the  sequeney  amplitude  threshold  Is  set  at  a different 
level.  A large  S/N  ratio  Improvement  results  la  each  case.  This,  however, 
could  have  been  accomplished  by  simply  low  pass  filtering  the  Walsh  spectrum 
since  the  low  order  (0  - A)  terns  dominate  in  the  reconstruction  of  the  pulse 
Data  2-3,  2-6,  and  2-7  - These  data  plots  have  the  same  input  data  but 
different  threshold  values  are  used.  No  S/N  ratio  Improvement  was  obtained 
in  2-3.  In  2-6  the  pulse  falls  deeper  into  the  noise  and  in  2-7,  the  pulse 
is  mo  longer  present  sfter  processing  has  taken  place.  For  this  input  data, 
the  pulse  lies  deep  enough  In  the  noise  so  that  the  higher  Walsh  terms  have 
amplitude  values  which  are  in  the  same  amplitude  range  as  the  low  order 
Walsh  terms.  The  high  order  Walsh  functions  have  many  sero  crossings  within 
the  width  of  this  one  pulse.  They  contribute  to  the  corruption  rather  than 
the  enhancement  of  the  flat-topped  pulse  characteristics.  By  high  order 
Walsh  functions,  X am  referring  to  the  functions  of  order  greater  than  32. 
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Sequencles  above  this  value  rarely,  if  ever,  contain  any  pulse  Information. 
Since  the  high  order  terms  contribute  in  this  example  as  much  as  do  the  low 
order  terms,  we  cannot  expect  an  S/N  improvement  by  just  amplitude  thresh- 
olding. In  data  2-7,  the  amplitude  threshold  is  set  too  high.  As  a result, 
many  of  the  low  order  terms  are  not  processed  and  the  pulse  is  lost. 

In  general,  amplitude  thresholding  of  the  Walsh  spectrum  will  not  give 
any  improvement  or  increased  likelihood  of  pulse  detection  at  low  3/N  ratios. 


z-z  ««a 


Data  2-4 


Data  2-5 


Data  2-6 


Data  2-7 


SECTION  V 


Data  3 Description 

Both  low  pass  flltetlng  and  amplitude  thresholding  of  the  Walsh  spectrum 
are  performed.  The  accuracy  of  pulse  reconstruction  la  dependent  upon  the 
number  of  terms  passed  by  the  low  pass  filter  and  to  the  sequency  threshold 
level  setting.  The  sequency  low  pass  filter  setting  Is  varied  from  20  to  30 
terms  which  are  passed  and  contribute  to  the  reconstructed  data  of  line  4. 

Notes  on  Data  3 

Data  3-6  and  3-7  have  identical  input  data,  but  3-6  does  not  amplitude 
threshold  the  sequency  terms.  3-7  has  a sequency  threshold  of  .3;  this 
results  in  smoother  line  4 data  and  avoids  the  false  pulse  Indication  of 

3-6,  line  5.  At  S/N  ratios  of  -5  db,  a pulse  can  still  be  filtered  from  the 

i 
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SECTION  VI 


CONCLUSIONS 

The  square  corner,  flat-copped  characteristics  of  Walsh  functions  make 
tb?m  Ideally  suited  for  the  processing  and  detection  of  rectangularly  shaped 
pulses  burled  In  noise.  Most  types  of  noise  are  sinusoidal  or  spiked  in 
shape.  If  a pulse  Is  buried  in  this  noise,  the  inverse  Walsh  transform 
requires  fewer  terms  to  reconstruct  a pulse  than  the  Fourier  transform  does. 

The  higher  the  order  of  a Walsh  function,  the  more  zero  crossings  It  will 
have.  Thus,  the  higher  erder  Walsh  aequencles  contribute  to  the  reconstruc- 
tion of  the  noise  rather  than  the  pulse.  By  low  pass  filtering  the  sequency 
spectrum,  much  of  the  noise  can  be  removed  while  Still  retaining  pulse 
Information.  Low  pass  spectrum  filtering  followed  by  amplitude  thresholding 
of  the  low  order  spectrum  elements  provlc.es  the  optimum  pulsed-slgnal 
detection  capability  at  low  S/N  ratios. 
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APPENDIX  I 


Software  Description 

All  of  the  Information  presented  In  the  data  sections  of  this  report 
was  obtained  from  computer  programs  written  In  Multics  Fortran.  These 
programs  were  Implemented  on  the  RADC  Multics  time  sharing  computer  system. 
The  FWT  subroutine  used  to  generate  the  Walsh  transform  can  be  found  in 
reference  3.  In  this  subroutine,  x is  the  256  point  array  containing  the 
real  data  samples  before  the  call  to  this  subroutine  is  made.  After  the 
subioutine  is  called,  x becomes  the  array  of  Walsh  coefficients.  Two 
successive  calls  to  this  subroutine  operating  on  x will  yield  the  original 
data  since  the  inverse  FWT  is  obtained  in  the  same  way  ac  the  forward  FWT. 

The  noise  used  in  the  generation  of  the  time  domain  input  data  was 
obtained  from  the  Multics  random  number  generator  through  the  call: 

call  random_$normal  (iseed,a) 

where  (a)  la  the  random  number  generated  and  lseed  is  the  starting  seed 
for  the  generation  of  the  random  noise.  This  noise  has  an  approximately 
Gaussian  distribution. 

The  graphics  were  produced  on  a Tektronix  4002  CRT  graphic  computer 
terminal.  Many  of  the  Tektronix  Plot  10  subroutines  were  used.  These 
subroutines  depend  on  the  output  of  certain  ASCII  control  characters  which 
the  Multics  system  would  not  output  directly  at  the  time  this  report  was 
written.  Thia  made  it  necessary  to  write  a PL/1  subroutine  utilizing  the 
Multics  users  subroutines  of  los_rawoutput  and  changemode.  This  subroutine 
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is  called  TOITTPT  and  it  will  output  all  control  characters.  The  Plot  10 
graphics  could  then  be  used  to  plot  the  data.  TOUTPT  does  not  pertain 
directly  to  this  report.  It  is  included  here  as  reference  material  for 
users  of  the  Multics  system. 

In  this  subroutine,  It  is  the  desired  ASCII  control  character  in 
decimal  number  form. 


toutpt:  proc(k)  ; 

del  ios_$write_ptr  entry  (ptr,fixed  bin, fixed  bln); 
del  buffO)  bit (36)  aligned; 
del  k bit (36); 

del  ios_$changemode  entry (char (*) ,char(*) ,char (*) ,blt(72)aligned) ; 

del  nlneblts  bit (9); 

ninebits-substr(k,28,9) ; 

del  status  bit(72)  aligned; 

del  gmode  char(32); 

del  sysprlnt  file; 

call  ios_Schangemode("u8er_i/oM,''rawo", gmode.  status) ; 
buff(l)-(ninebits)  ||  ”000000000"b  ||  "000000000"b  1 1 "000000000,,b ; 
call  ios_Swrite_ptr  (addr(buf f (1)) ,0,1) ; ” 

cull  ios_$changemode("user_i/o", gmode, gmode, status) ; 
end  toutpt; 


i 
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MZTMC  8YSTD4 


IASI  UNITS: 

Quantity  _ Unit 


length 

metre 

mate 

kilogram 

lima 

tacond 

electric  currant 

ampere 

thermodynamic  temperature 

kalvin 

amount  o f aubatanca 

mole 

lumlnoua  tnteneity 

candela 

SUmJMXNTAXT  UNITS: 

plana  an(la 

radian 

aolid  angle 

at  arad  Ian 

DUUVni  UNITS: 

Acceleration 

metre  par  mcond  equated 

activity  (of  a mdioactlva  aouica) 

dlilntagratlon  par  aacond 

angular  accalaratton 

radian  par  aacond  equated 

angular  valocity 

radian  par  aacond 

araa 

aquara  mitre 

danaity 

kilogram  par  cubic  matte 

electric  capadtanca 

(wad 

electrical  conductance 

tlemena 

electric  fiald  atrangth 

volt  par  matra 

eiectrk  tnductanca 

henry 

alactrtc  potential  difference 

rail 

alactrtc  raaiatanca 

ohm 

alactromottva  (area 

veil 

anargy 

Me 

antrapy 

|oule  pm  kalvin 

force 

newton 

frequency 

hard 

llluminanca 

hn 

luminanca 

candela  par  aquara  metre 

lumlnoua  flua 

lumen 

magnatlc  field  atrangth 

ampara  par  metre 

magnetic  flun 

webar 

magnatlc  flun  danaity 

tea  la 

magnetomotive  (area 

ampara 

po  war 

watt 

praaaura 

paacal 

quantity  of  alactilclty 

coulomb 

quantity  of  hunt 

|oula 

radiant  intanaity 

wort  par  Wared  Ian 

•pacific  haat 

ioule  par  kHograra-katvIn 

atraaa 

paacal 

tharmal  conductivity 

watt  par  metrv-kelvla 

velocity 

■MtW  ptf  Mcond 

aiacoalty.  dynamic 

paaral-aacond 

viacaaMy.  kinematic 

aquara  metre  par  aacond 

voltaga 

vail 

volume 

cubic  matra 

wavenumber 

reciprocal  matra 

work 

|aum 

stpuraui 

_ Multiplication  factor* 

1000  000  000  000 

- 10*» 

1000  000  000 

- 10* 

t 000  000 

* 10* 

1000 

- 10* 

100 

• 10* 

10 

- 10' 

0.1 

- 10-' 

001 

« 10- • 

0 001 

• 10-' 

II  OIK)  001 

-10  • 

o.oou  oco  ooi 

- io-» 

o.ooo  ooo  ooo  ooi 

• to-" 

0 000  000  000  coo  001 

- 10-" 

0.000  000  000  000  000  001 

- I0-" 

* To  ha  avoided  where  pnaalbla 

